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ABSTRACT 

One hundred and forty six long-period red variable stars in the Large Magellanic Cloud 
(LMC) from the three year MOA project database were analysed. A careful periodic 
analysis was performed on these stars and a catalogue of their magnitudes, colours, 
periods and amplitudes is presented. We convert our blue and red magnitudes to K 
band values using 19 oxygen-rich stars. A group of red short-period stars separated 
from the Mira sequence has been found on a (logP, K) diagram. They are located at 
the short period side of the Mira sequence consistent with the work of Wood and Sebo 
(1996). There are two interpretations for such stars; a difference in pulsation mode or 
a difference in chemical composition. We investigated the properties of these stars to- 
gether with their colour, amplitude and periodicity. We conclude that they have small 
amplitudes and less regular variability. They are likely to be higher mode pulsators. 
A large scatter has been also found on the long period side of the (logP, K) diagram. 
This is possibly a systematic spread given that the blue band of our photometric 
system covers both standard B and V bands and affects carbon-rich stars. 

Key words: catalogues - stars: AGB and post-AGB - stars: oscillations - stars: 
variables:other - Magellanic Clouds 
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1 INTRODUCTION 

There are many types of variable stars whose mechanism for 
variability (mainly pulsation, rotation, eruptive processes, 
and geometric effects) is well studied and well understood 
observationally and theoretically. Among them, the pulsat- 
ing variable stars have been investigated vigorously because 
they yield information on their stellar mass and their inte- 
rior physical properties. They are a powerful tool for mea- 
suring distances in the Galaxy and extra-galactic systems 
(1) because they are luminous and (2) because they have 
been shown empirically to obey a period-luminosity (PL) 
relation. 

The PL relation for classical Cepheids was discovered 
by Leavitt (1908) and used to measure the distance of galax- 
ies. However, there is a discrepancy between the distance to 
the LMC derived from the PL relations for Cepheids and 
that derived for RR Lyraes (Walker 1992) . This discrepancy 
may be due to uncertainties in the absolute calibration and 
a sensitivity to metallicity effects (Caputo 1997). It is im- 
portant to establish PL relations of any new objects which 
are independent of these variable stars. 

A PL relation for long period red variables (Miras and 
semi- regular variables (hereafter SRs)) will be a useful dis- 
tance indicator because these stars are numerous and have 
high luminosities. However, they are the least studied star 
type because of their complex pulsation mechanism. Diffi- 
culties arise from uncertainties in our understanding of the 
convective energy transport and the opacity resulting from 
the many molecular absorption lines. 

We present in this paper the results of our study of LMC 
Miras and SRs using the MOA database. This database was 
generated from the MOA project microlensing survey. 

A PL relation for Miras of the LMC was presented by 
Glass & Lloyd-Evans (1981) and Glass & Feast (1982). They 
observed 11 Miras in the infrared bands J,H,K. However, 
there was a large amount of scatter in the plot {a ~ 0.25 
mag) because their observation time was limited and it did 
not cover the whole phase of variations. 

To reduce the scattering in the plot due to poor phase 
coverage, extensive infrared observations have been made 
(Feast et al. 1989; Hughes & Wood 1990), and a narrow PL 
relation was obtained (cr ~ 0.13 mag). Feast (1996) showed 
that the Miras in the LMC follow a narrow PL relation in 
both absolute bolometric magnitude, Mi,oi , and absolute K- 
magnitude, Mk, and he concluded that all Miras pulsate in 
the same mode. 

Wood & Sebo (1996) found two sequences in a, K band 
PL diagram for the small amplitude SRs in the LMC. One 
was the well-known Mira sequence, and the other (eight 
stars) was located at the position which is parallel to the 
Mira sequence and their periods were approximately half 
the periods of Miras on the PL diagram. This discovery 
encouraged a hypothesis that the pulsation mode of these 
stars in the secondary sequence must be in a higher mode. 
Using the Hipparcos database of nearby stars. Bedding & 
Zijlstra (1998) discovered the SRs which are located in the 
secondary sequence of Wood & Sebo (1996). Bedding and 
Zijlstra interpreted the sequence of SRs as evolutionary (see 
Whitelock (1986)). 

Feast (1999) suggested that both Miras and SRs pul- 
sate in the same mode, possibly both in the first overtone. 



The arguments were based on the angular diameter mea- 
surements of nearby Miras and the PL relationship of the 
Miras and SRs in globular clusters. He agreed with Bedding 
and Zijlstra (1998) that the sequence which is located on 
the shorter period side of the Mira sequence in the PL plane 
is an evolutionary track ending at the Mira sequence. Feast 
also suggested that the Mira sequence would be a collection 
of variable stars with different chemical abundances and/or 
with different stellar mass. In this case, it was not necessary 
that Miras and SRs pulsate in different modes. 

As described above, there are two theories on the short- 
period stars in the PL diagram for long-period red variable 
stars. Considering that the pulsation features of the Miras 
and SRs depend on several factors, other physical properties 
must also be considered. The study of stellar colour will be 
interesting because this gives us the temperature of stellar 
surface, stellar radius and mass which are the important 
factors of the pulsation period. 

Recently, the continuous observation of a great number 
of stars using CCD cameras which cover a large image area 
became possible. The recent report of Wood (2000) based 
on the K band photometry and the period derived from the 
MACHO database indicated plausible sequences on the PL 
diagram. This diagram shows at least five sequences falling 
parallel to the Mira sequence. According to Wood, the three 
sequences containing the Mira sequence were explained as 
being different pulsation modes. However, more information 
is required to settle this issue. The features of the other two 
sequences are unknown, and it is expected that further inde- 
pendent observations and detailed investigation will reveal 
the nature of these stars. 

We worked on the problem using the MOA database 
which contained data for a large number of stars in the LMC. 
The period, magnitude colour, amplitude and periodicity for 
each star was determined. 

We found an independent group of short period stars 
in a sample of 146 stars. We investigated the properties of 
these stars through their colour, amplitude and periodicity, 
and discuss the possibility of that they are higher mode pul- 
sators. 

Following an introduction to the MOA project, we de- 
scribe our observation and photometry systems in Sections 
2 and 3. To obtain the light curves of variable stars, a care- 
ful calibration of each image frame is required. This pro- 
cedure is described in Section 4. In Sections 5 and 6, the 
data reduction and period analysis methods are described. 
In Section 7, the colour-magnitude diagram (CMD) of the 
LMC stars observed using our system is presented and the 
period-magnitude relation will be given in Section 8. The 
period-colour relation is presented and discussed in Section 
9. The result of the conversion of MOA magnitudes to K 
band values will be given in Section 10, along with a K, log P 
diagram. In Section 11 we discuss the possibility of the ex- 
istence of multiple sequences in the (log P, K) plane with 
the amplitude distribution. In Section 12 we summarise our 
results. 



2 THE MOA PROJECT OBSERVATIONS 
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Figure 1. The transmission of the filters and quantum efficiency 
of MOA-cam2 are shown as a function of the wavelength. We 
used two broad-band filters {MOAb, MOAh), where MOAg 
covers standard B and V, and MOAji covers the Kron-Cousins 
RI region. 

2.1 Observational Instruments 

The MOA project commenced observations in May 1996 
with the aim of observing microlensing events towards the 
LMC, SMC and Galactic Bulge. The observational program 
was undertaken at the University of Canterbury's Mt. John 
Observatory (44° 5", 1030m above sea level) in the centre of 
the South Island of New Zealand. 

The MOA project uses a 61cm Ritchey- Chretien 
Cassegrain telescope with modified optics to give a wide 
FOV at f/6.25. 

For the first set of observations in 1996 (Series 1), we 
used the original f/13.5 telescope optics and monitored ~ 
3 X 10^ stars in three LMC fields with the mosaic CCD 
camera, MOA-caml. MOA-caml was constructed with nine 
TI TC215 CCDs with 1000 x 1018 pixels whose pixel size was 
12^m X 12fim. From January 1997 to July 1998, in Series 2, 
the optics were altered to f/6.25. MOA-caml covered 3 fields 
in the LMC bar (1 x 10'^ stars) and 2 fields in the SMC 
(~ 4 X 10^ stars). The current set of observations. Series 
3, started in August 1998. For this series, a new camera, 
MOA-cam2 is in use. M0A-cam2 consists of three abutted 
SITe 2047 x 4095 pixel thinned, back-illuminated CCDs with 
15/im xl5/im pixels. The pixel size is 0.81", giving a field 
of view of 0.5° x 0.9° for each chip. The details of MOA- 
cam2 are in Yanagisawa et al. (2000). In this series we are 
observing 16 LMC fields (4.4 x 10^ stars), 8 SMC fields (9.3 x 
10^ stars) and 14 Galactic Bulge fields. 

We used two broad-band filters to maximise photon col- 
lection. The MOA blue fifier (MOAb) covers 395 - 620nm 
(Aeff ~ 500 nm) and MOA red fiher (MOAr) covers 
620 — 1050 nm (Aeff ~ 700 nm). The spectral transmission 
curves for each filter are shown in Fig. 1. 

The MOA instrumental photometry has been trans- 
formed to the standard Johnson V and Kron-Cousins R 
(Rkc) by comparisons with about one hundred stars im- 
aged by Hubble Space Telescope (HST). We describe these 
details in Section 4. 



The variable star study based on the MOA database 
was reported in Hearnshaw, Bond, Rattenbury, et al. (2000) , 
Noda, Takeuti, Bond, et al. (2000), and Takeuti, Noda, 
Bond, et al. (2000). 



2.2 OBSERVATION DETAILS 

Regular 300 second exposures were made towards the Large 
Magellanic Cloud and the Small Magellanic Cloud. 180 sec- 
ond exposures were made for the Galactic Bulge (GB) fields. 
Under the current system, it takes 70 seconds to read out 
the 3 CCD images for each exposure. Dark frames for each 
exposure were prepared every photometric night and dome 
flats or dark-sky fiats were also frequently taken. 

Observational data were archived on DLT tapes and all 
the data were analysed in both Japan and New Zealand. 
Series 1, 2 and 3 images were analysed using DoPHOT pho- 
tometry software (Schechter et al. 1993). The size of the 
Series 2 database created by the subsequent photometry is 
~ 40 Gbytes. 

A new online analysis system based on image subtrac- 
tion has been developed in New Zealand. In the year 2000 
this system was able to find several microlensing events in 
real time and issue alerts to the astronomical community 
via the WWWQ. We are re-analyzing some parts of our data 
with this method. 



3 DATA ANALYSIS 

Before undertaking data analysis using the DoPHOT pro- 
gram, all image files were prepared with dark and the flat 
field frames. DoPHOT was used in a fixed position warm 
start mode, using the best photometric images as templates. 
It takes about 50 minutes to analyze one LMC frame with 
our computer. There were more than 400 measurements in 
the Series 2 (1.6 year) and Series 3 (1.4 year) database. We 
obtain only relative magnitudes from the DoPHOT analysis. 
We also use non-standard filters, thus a transformation of 
our instrumental magnitude to the standard V and Rkc is 
required. This transformation is described in the next sec- 
tion. 



4 CALIBRATION FOR MOA PHOTOMETRY 

The transformation of the MOA photometry to the stan- 
dard system was carried out using Hubble Space Telescope 
(HST) image data. More details are given by Kato (2000). 
We defined the transform relation for MOA-caml and MOA- 
cam2 separately because their CCD quantum efficiencies 
were quite different (the quantum efficiency of MOA-cam2 
was twice as high as that for MOA-caml). Although the fil- 
ter system is exactly the same, this will cause serious differ- 
ences in the relative magnitudes returned by DoPHOT. The 
DoPHOT magnitudes were calculated using photon quan- 
tities only, so a difference in CCD quantum efficiency will 



t http:/ /www. phys.canterbury.ac.nz/^physib/alert/alert. html 
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affect the DoPHOT output directly. To get standard magni- 
tudes for MOA-cam2, we performed the procedure described 
below. 

Photons coming from a star show a complicated scatter- 
ing pattern around a central position caused by atmospheric 
turbulence. A point-spread function (PSF) models the scat- 
ter. The DoPHOT code estimated each magnitude using an 
elliptical Gaussian with seven free parameters as a model 
for a stellar PSF. A common PSF was adopted for all stars 
(bright stars and faint stars) in a frame, however it was clear 
that there were significant differences between the profiles 
for bright stars and that for faint stars. The former were 
weighted by the photon statistics of the star itself, and the 
latter were weighted by the photon statistics of the sky, so 
the latter should show larger scattering. The detailed study 
of the PSF indicated that the function often differs from the 
elliptical Gaussian ( Alard & Lupton 1998) . The error due to 
this discrepancy is approximately ±0.05 mag for the bright 
stars and ±0.09 mag for the faint stars. It was not a criti- 
cal problem if we define a mean magnitude with more than 
100 measurement points to study stellar intensity, because 
the scatter due to PSF model inadequacy usually occurs at 
random. 

As we were using large CCD chips, we have to deter- 
mine the position dependent error for each CCD chip. In 
order to investigate this, we observed the globular cluster 
NGC3201 at various positions on each chip, and compared 
stellar magnitudes with published magnitudes obtained us- 
ing B,V,R,I filters by Alcaino et al. (1989). We found that 
the error was smaller (±0.015 mag) than the DoPHOT PSF 
photometry error (< ±0.09 mag). This was a systematic er- 
ror however, and did not disappear even if we deal with a 
mean magnitude. 

In the third step, also using the NGC3201 data, we 
computed the chip to chip differences for the magnitude zero 
point: 



MOAB(^chip3) = MOAB(^chipi) - 0.23 
MOAr^mpS) = MOAn(Mpi) - 0.31 
MOABiMp3) = MOAb(cMp2) - 0.34 
MOAr^mpS) = MOAr^mp2) - 0.41 



(1) 
(2) 
(3) 
(4) 



We used the chip3 scale as our standard as it gave the 
most linear response. 

For the final step, we converted our MOAb and MOAr 
magnitudes to standard Johnson V and Kron-Cousins R 
(Rkc) respectively. The brightness of stars was acquired 
from an ADU (Analog Digital Unit, 1 ADU corresponded 
to 1.3 electrons for M0A-cam2) value which corresponded 
to the number of photons measured with the assumed PSF 
via the DoPHOT photometry. We define the zero-point of 
magnitude for MOA as 1 ADU. Therefore all of the ob- 
tained magnitudes are negative, and we have to correct the 
zero-point of our magnitude in order to obtain actual mag- 
nitudes. We analyzed the HST data of the LMC obtained 
with the Wide Field Planetary Camera2 (WFPC2) to cor- 
rect for atmospheric seeing. Since the filters for WFPC2 had 
been transformed to the Johnson V and Rkc by Holtzman 
et al. (1995), we were able to compare the converted HST 
magnitudes directly with those of our one hundred stars. 
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Figure 2. Colour-Magnitude Diagram (CMD) for the 146 stars in 
our sample with 9,050 randomly chosen stars (dots). The vertical 
axis indicates Rm- All 146 stars are located on the tip of Giant 
Branch, which means they are in the AGB stage. The seven short 
period stars are indicated with + symbols. The ffi mark indicates 
nlmc3-214266 (see Section 7). The other stars are plotted as open 
circles. 



Finally, we found a relation to convert the Vm and Rm 
for MOA-cam2; 

Vm = MOAB-0.16{±0m)*{MOAB-MOAR)+27.41{±0m){5) 

Rm = MOAr+0.29{±0.01)*{MOAb-MOAr)+27.24{±0.01){6) 

Colour{Vm-Rm) = Q.^2,*Colour{MOAB-MOAR)+0.l?,{7) 

For MOA-caml; 

Vm = MOAs-0.10(±0.01)*(A/O4s-MOAb)±25.55(±0.01)(8) 

Rm = MOyls±0.31(±0.02)*(A/OAB-AfOAfl)+25.45(±0.1)(9) 

Colour(Vm-Rm) = 0.59*Co;our(A'/O^s-A'/O^fl)±0.10(10) 

Hereafter, we deal with Vm and Rm as the MOA stan- 
dard magnitudes. In Fig. 2, we show the Colour Magnitude 
Diagram (CMD) obtained by MOA-cam2 for our observa- 
tion fields in the LMC bar. The limiting magnitude was 
around 20.5 mag for MOA-cam2, and 19 mag for MOA- 
caml. Since our passbands were broad and covered frequen- 
cies from violet to near infra-red bands, the colour is not 
sensitive for blue stars. 
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To minimize the effects of the non-uniformity of the sky 
and system on the raw magnitudes (i.e. DoPHOT output 
magnitudes) we: 

(i) Divide the CCD chips into smaU regions (200 x 200 
pixels) 

(ii) Select eighty bright stars which had more than 100 
measurement points in each small region 

(iii) Regard stars with a small scatter as 'non- variable' 
stars, i.e., where the number of observed points more than 
4 sigma from the average magnitude was less than 15 

(iv) Calculate magnitudes for all stars by comparison 
with the 'non-variable' stars. That is, the MOA standard 
magnitudes were determined relative to the magnitudes of 
the near 'non-variable' stars in each small region. 

Through the revision of our calibration procedure, we 
obtained progressive and interesting diagrams compared 
with our previous reports (Hearnshaw, Bond, Rattenbury, 
et al. 2000; Noda, Takeuti, Bond, et al. 2000; Takeuti, Noda, 
Bond, et al. 2000). Details are given in later sections. 



5 DATA REDUCTION 

We selected long-period red variables from the Series 2 
database using the following criteria: 

(i) Each star should be 'photometric' and have a 'single 
star shape' as judged by DoPHOT. 

(ii) Light curves should be sufficiently long (> 100 
points) . 

(iii) Light curves of each star must have large variability, 
i.e., the reduced for a fit to mean brightness should be 
larger than Xvm > Xo and xlm > Xo- Here > Xo corre- 
sponded to the upper 10% limit of the reduced distribu- 
tion. For the MOA database Xo = 4.3. 

To obtain large amplitude variables, we applied the addi- 
tional criteria: 



(iv) Low amplitude stars at an early stage of the study 
are rejected by removing those stars which had a magnitude 
difference between the fourth highest and fourth lowest mag- 
nitude values of less than 1.3 magnitudes. This difference 
was larger than the amplitude SRm determined in Section 
6. This selection was done for the results of MOA-caml sys- 
tem without using the colour dependence in equations (8) - 
(10). This sample of large amplitude stars might be affected 
when we use the correct colour relation. 

(v) We select redder stars for large amplitude; {Vm — 
Rm) > 0.4 for the MOA-caml system. 

We obtained 1,067 light curves for the candidates of 
variable stars with large amplitude after applying criteria 
(iv) and (v) on the Series 2 database. Light curves from Se- 
ries 2 and Series 3 were cross referenced to stars from their 
respective templates. We rejected unreliable light curves, 
and finally obtained 908 light curves. In Fig. 3 we show an 
example of a light curve. 

We also found 485 small amplitude variable star can- 
didates by applying a further criterion after numbers (i) to 
(iii) above. We analyse the large amplitude variables in this 
work. We will present the small amplitude variable star re- 
sults in a future paper. 



14 


. 1 1 1 1 1 1 


nlmol-280227- 


-15 






16 






17 


I.I.I 


1 . 1 



400 



600 



BOO 



1000 



1200 



1400 



15 
«16 



17 



' 1 ' 1 ' 1 


1 1 1 1 1 . 

(a) nlmol-324832 : 






r . 1 . 1 . 1 


1 . 1 



400 



600 



800 



1000 



1200 



1400 



nlmc3-B8048 



_i_ 



_i_ 



_i_ 



_i_ 



_i_ 



14 
3 15 
-16 
17 

400 600 800 1000 1200 1400 

JD-2450000 

Figure 3. Examples of light curves for Rm - The periods for these 
stars were identified as 268.0, 133.5, and 131.2 days respectively 
by the PDMM and PDM algorithms (see Section 6). The time 
axis spans from January 1997 to December 1999. The error of 
each observation point is the DoPHOT PSF error. The stars 
nImcl-324832 and nlmc3-88048 were a short-period red star and 
a short-period blue star discussed in Section 11. Lack of data af- 
ter JD = 2451400 for nlmcl-280227 and nlmcl-324832 were due 
to instrumental problems with chip 1 on MOA-cam2. 

6 PERIOD ANALYSIS (PDMM) 

As the first stage in the period analysis, we used the Phase 
Dispersion Minimization method for MOA (PDMM) which 
is a PDM method (Stellingwerf 1978) modified to operate 
on large photometry datasets. The PDMM algorithm finds 
a likely light curve period together with a reliability param- 
eter, 0, for a large number of light curves automatically. 
9 — indicates a complete cyclic change and 9 = 1 indi- 
cates no periodicity. The amplitude, the difference between 
the magnitudes of the brightest and the faintest bins, was 
also calculated. We converted the magnitudes of each bin 
to intensities and calculated the mean intensity. We denote 
the mean magnitude values as (Vm) and (Rm)- We executed 
this code using both colours {Vm,Rm) for the 908 selected 
light curves. 

Since our dataset included scattered data points, the 
PDMM algorithm sometimes gave integral multiples or sub- 
multiples of the period, or in the worst case an incorrect 
result. To avoid simple errors in the periods obtained by the 
automatic analysis, we checked all code output carefully, and 
ran the PDM code manually to determine the most likely 
periods. We accepted a period which had good agreement 
between both colours. In Fig. 4, we present an example of 
PDMM figures. The upper diagram shows the reliability of 
the testing period. A small value of 9 indicates a likely value 
for the period. The bottom diagram shows the light curve 
folded with the most likely period. 

The total dataset (Series 2 to Series 3) extends over 
about 3 years (~ 1,060 days), so periods less than ~400 days 
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Figure 4. A sample of PDMM algorithm output for the same 
star, nlmcl-280227, as shown in Fig. 3. The top panel shows the 
reliability of period determination. The 9 is the sum of the dis- 
persion of data in each bin for the given period. The period corre- 
sponding to the minimum of 8 is chosen as the most likely period. 
The bottom panel shows the folded phase curve with the most 
likely period. The error bars indicate the standard deviation in 
each bin. The period of this star was determined as P = 268.0 
days with 6 = 0.21. 



should be well determined because these have more than 
2.5 cycles of variation within the observation period. We 
classified 908 candidate hght curves by their features of pe- 
riodicity. Among 908 stars, 329 stars show clear variability 
and of these, 259 stars had periods which were relatively 
well-determined. We rejected hght curves like eclipsing bi- 
naries, then chose hght curves whose periods were longer 
than logP > 1.5 (30 days). We also required logP < 2.6 
(400 days) reflecting the dataset time span. Finally, 146 light 
curves were left. The reliability parameter 6 for these curves 
is less than 0.67 (the mean ~ 0.38) with 15 folding bins. We 
show the distribution of 9 in Fig. 5. 

Other candidates, the stars with poor determined pe- 
riod, short period or small amplitudes, will be the subject 
of future work. Here we give observational values for 146 
stars in Table 1 (ID, a, 5, (Rm), (Vm), SRm, SVm, Period, 
6). ID is the MOA database star number, (a, 5) is the right 
ascension and declination respectively. Comparing the co- 
ordinates of 629 stars in Tables 2 ~ 6 of Hughes & Wood 
(1990) with our sample, we identified at least 34 stars as a 
control. The ID numbers of their tables are indicated in the 
Remarks column. Spectral types given in their tables are 
also shown. The 34 stars consisted of 21 oxygen-rich stars 
and 13 carbon-rich stars. In Fig. 6, we show the relation 
between the period determined via the PDM method and 
the period tabulated in Hughes & Wood (1990) for 30 stars 
whose periods are well-determined by us. It is clear that the 



53 O 



146 stars 



0.2 



0.4 



0.6 



0.8 



Figure 5. The distribution of 9 calculated with 15 bins in 
PDMM. The mean of 9 is approximately 0.38 for the 146 samples. 
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Figure 6. The comparison between the periods determined by 
our PDM (horizontal axis) and the period tabulated in Hughes 
& Wood (1990) (vertical axis). The differences are quite small 
(< 2%). 



periods are coincident and the differences between MOA and 
Hughes & Wood (1990) is less than 2 %. We used these 30 
stars for the conversion to the K magnitude system as de- 
scribed in the following section. Here, we note that Hughes 
& Wood (1990) tabulated two periods, 170 days and 319 
days, for the star nlmc3-288533. We adopted 319 days as 
their period because we obtained 307.6 days for the period 
of this star via our analysis using the PDM method. 
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Table 1. 146 selected long-period red variables 
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Table 1 - continued 146 selected long-period red variables 
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Note: (Rm), {Vm) are the intensity mean magnitude of Rm or Vm respectively, and 9 is the PDM algorithm reliability parameter. 
(a)~(g) and (1) are the stars referred to in the text. Numbers in the Remarks column are the identification number (with spectral 

type) in Hughes and Wood (1990). 



7 COLOUR MAGNITUDE DIAGRAM (CMD) 

Fig. 2 shows the distribution of our sample of 146 variable 
stars on a colour-magnitude diagram (CMD) with 9,050 ran- 
domly chosen stars for reference. In this figure, the 146 stars 
were located at the tip of the Giant Branch. According to 
the stellar evolution theory, stars with initial mass of about 
1 — 8Mq evolve into Asymptotic Giant Branch (AGB) stars 
with carbon and oxygen cores. Although we only required a 
large variation in each light curve and restricted the colour 
condition and period in our selection criteria for variable 
stars (Section 6), these stars were located in one particular 
region in the CMD. Consequently we presumed that these 
stars form a group which were in a common evolutionary 
stage and pulsate by a common mechanism. However there 
was an extremely bright star (indicated by a © mark) in this 
CMD, nlmc3-214266, which we have mentioned below. Also 
we will describe the seven stars with -I- marks in Section 11. 
nlmc3-214266 and these seven stars are denoted by (1) or 
(a) ~ (g) respectively in Table 1. 



8 PERIOD-LUMINOSITY RELATION 

We show a plot of (Vm) magnitude versus period in Fig. 7. 
Open circles indicate stars whose Rm amplitudes (SRm) are 
smaller than 2.2 mag. Filled circles indicate the stars whose 
amplitude are larger than 2.2 mag. The line indicates the 
regression fit for the stars with log P > 2.2; 

{Vm) = 2.01(±0.43) * log P + 12.10(±1.05) (11) 



o amplitude(Rm) < 2,2 mag 

• ainplitude(Rm) > 2.2 mag 

ffi nlmc3-214266 

+ short period stars (7 stars) 




logP[days] 

Figure 7. The Period-Magnitude {(Vm)) diagram. Open circles 
are the stars whose amplitude {SRm) are less than 2.2 mag. Filled 
circles are the stars whose amplitude are larger than 2.2 mag. The 
extremely bright star, nlmc3-214266 is indicated by the © mark. 
This star is much brighter than the other stars with similar period 
in this diagram. The cross marks indicate the stars of group B 
(indicated in Section 11). The downward line indicates the least 
square fitting for the stars with logP > 2.2. 



The trend to be fainter in {Vm) has been known in pre- 
vious studies. The star, which was extremely bright in the 
CMD, nlmc3-214266 {{Vm) = 14.69, logP = 2.05, © mark). 
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Figure 8. The Period-Magnitude {(Rm)) diagram. See Figure 7 
for an explanation of the symbols used. 



was also much brighter than the other stars with similar pe- 
riods in this diagram. We also show the (logP, (Rm)) plot 
in Fig. 8. The slope of this distribution for the stars with 
log P > 2.2 was less than that in {{Vm) , log P) plane but still 
had a downward tendency. The best fit line for these stars 



(Rm) = 1.13(±0.33) * log P -I- 12.84(±0.80) 



(12) 



The star, nlmc3-214266, was again quite bright {{Rm) ~ 
13.80) and was separate from the other samples. The scatter 
of {Vm) magnitude {a ~ 0.52) was larger than that of {Rm) 
{a ~ 0.39), which was possibly due to the difference between 
oxygen-rich stars and carbon-rich stars. We will discuss this 
in Section 10. 



9 PERIOD-COLOUR RELATION 

We show the period-colour {log P,{Vm) — {Rm)) plot in 
Fig. 9. The upward line indicates the regression line for the 
stars with log P > 2.2; 



{Vm) - {Rm) = 0.88(±0.16) * log P - 0.74(±0.38) 



(13) 



In this figure, the colour of most stars with longer pe- 
riods (log P > 2.2) were much redder than that of the short 
period stars, and the relationship was more constrained 
((5 ~ 0.19) compared with that shown in Fig. 7 and Fig. 8. 



9.1 nlmc3-214266 

In Fig. 9, star nlmc3-214266 is positioned in the main group. 
Comparing with the other stars of similar period (open and 
filled circles), the effective temperature of this star should be 
similar since the colour {{Vm) — {Rm)) was almost equivalent. 
The brightness of nlmc3-214266 was high in both {Vm) and 
{Rm), therefore its radius must have been large compared 
with the stars in the main group. The star was almost cer- 
tainly HV 1011 which was a suspected Cepheid (Kurochkin 
1992). Kurochkin was unable to obtain a period for this star 
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• ainplitude{Rin) > 2,2 mag 
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Figure 9. The Period-Colour {(Vm) — {Rm)) diagram. The rela- 
tionship shows a narrower band compared with Fig. 7 and 8. The 
seven stars with short period (+ symbols, see Section 11.2) show 
the extreme red colour, while nlmc3-214266 is not conspicuous 
(see Section 9.1). 



and our results do not conclusively suggest that this star is a 
classical Cepheid. When the pulsation period is the same, a 
larger star must have a higher mass considering the relation 
between the period and the mean density of stars. Therefore 
we assumed that nlmc3-214266 was a massive star evolved 
from the main sequence, and not an AGB star. The light 
curve is presented in Fig. 10. The star nlmc3-214266 was a 
luminous semi-regular variable. In the Galaxy, V810 Cen- 
tauri is one of the variable stars which is more luminous 
and has a higher temperature than the AGB in the CMD. 
Kienzle, Burki, Burnet, and Meynet (1998) estimated that 
the spectral type of this star was F8 la, and the period was 
approximately 156 days, Mboi = —8.5, and the initial mass 
was ~ 25 Mq. Such a heavy star will become a supernova 
rather than a AGB star (e.g. Becker, 1998). 

9.2 Short Period Red Stars 

In Fig. 9, several stars whose periods were log(P/ci) < 2.15 
were distributed towards the red region and were separated 
from the main group of stars. If the colours are different 
between stars with a similar period and luminosity, their 
mass or pulsation property {Q = Py/p) should be different. 
It is quite important to investigate their luminosity because 
the relation between luminosity and colour gives informa- 
tion about their radius, i.e., mass and pulsation parameter 
(Q) of each star. Although our {Vm) and {Rm) are trustwor- 
thy they cover too small a wavelength range to enable the 
energy distribution of such red stars to be discussed directly. 
Therefore we study the relationship between our magnitudes 
and the near-IR luminosity in Section 10. 



9.3 nlmcl-239264 

We note that there were no data points for the minima of 
the light curve for nlmcl-239264. Therefore the mean magni- 
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Figure 10. The light curve of nlmc3-214266. This star is quite 
bright in the period-magnitude diagrams (Fig. 7, 8), not red 
in the period-colour diagram (Fig. 9), and again bright in the 
(logP, i^m) diagram (Fig. 15). 



tudes, (Vm) and (Rm), for this star were calculated by using 
only the bright phases. This gives just the brighter limit of 
the mean magnitude for this star. Regarding the period, we 
assumed that the determined period for the star (P — 129.5 
d, or logP = 2.11) was accurate because there were three 
maxima in the light curve. On Figs. 7 and 8, the star was 
located near the extension of the least squares fitted line for 
stars with log P > 2.2. When we used a dataset covering 
complete cycles, the stars were located at fainter positions. 
We used (Vm) — {Rm} ~ 1.44 as the blue or lower limit of 
the colour average in Fig. 9. The star was located about 
half a magnitude redder than the least squares fitted line 
in Fig. 9. Considering the incompleteness of our dataset, we 
may assume that the star nlmcl-239264 should be located 
at a fainter position on Figs. 7 and 8, and at a redder or 
higher position on Fig. 9. 



9.4 nlmc3-306907 

The star, nlmc3-306907, is the remarkable star because of 
its unique colour {{Vm) — {Rm)) variation, i.e., the colour 
was bluer when {Vm) and {Rm) were fainter. The light vari- 
ation is shown in Fig. 11. Among 146 samples, nlmc2-151839 
and nlmc3-191203 also showed the same property, though to 
a lesser extent. The colour variations of all the other stars 
were fainter when the (Vm) and {Rm) were brighter, which 
is a common feature of radial pulsators. Such an extraordi- 
nary property might indicate the development of a particular 
molecular band caused by a different chemical composition 
in the stellar surface. 
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Figure 11. The light curve of nlmc3-306907. The colour {{Vm) - 
{Rm)) is bluer when the {Rm) and {Vm) are fainter. 



10 PERIOD TO A'-MAGNITUDE RELATION 

It is well known that Miras fall on a narrow PL sequence 
on the (log P, K) plane (Feast et al. 1989; Hughes & Wood 
1990). Wood & Sebo (1996) discovered a second sequence 
located on the short period side of Miras, while Bedding 
& Zijlstra (1998) found strong evidence that Semi-Regulars 
(SRs) near the Sun's position in the {\og P,K) plane lie in 
a region which coincides with the second sequence found by 
Wood & Sebo (1996). Also the existence of at least two (or 
four) sequences other than Miras was indicated from the 
analysis of MACHO database (Wood 2000). 

The iC'-band (~ 2.2/im) corresponds closely to the max- 
imum radiation fiux of long-period red variables, so that the 
discrepancy between /i"- magnitudes and bolometric magni- 
tudes should be relatively small. Because previous studies 
have been done with K magnitudes, we have to convert our 
observational magnitude {Vm),{Rm) into Tf- magnitudes in 
order to compare our results with the established PL rela- 
tion. 

First, to define our conversion scale, we identified 34 
stars which were the same as those published by Hughes & 
Wood (1990) from our 146 sampled long-period red variable 
stars (see Section 6). The 34 identified stars consisted of 
21 oxygen-rich (M-type) stars and 13 carbon-rich (C-type) 
stars. Furthermore, we selected 19 oxygen-rich stars whose 
periods were well-determined in our analysis, and used them 
to calibrate our Vm , Rm magnitudes to K values; 



{Rm) - Kshv = 2.10 * {{Vm) - {Rm)) + 1.22 



(14) 



Kshv was the iC-magnitude tabulated in Hughes & 
Wood (1990) {shv is the prefix of the star number). Here- 
after we write Km in place of Kshv as our custom K magni- 
tude. We discuss here a few reasons for the scatter (a ~ 0.15 
mag) in this relation. Firstly, and possibly most importantly, 
the if- magnitudes obtained by Hughes & Wood (1990) were 
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determined with only a single observation, whilst we were 
dealing with the average magnitudes of at least one hun- 
dred measured points. Although amplitudes of long-period 
red variable stars are large (~ 0.6 mag) in the J^-band, 
a single epoch observation provides the magnitudes on the 
various phase of light variations, so it should have caused 
extra scattering in our Km conversion. Secondly, our trans- 
formed magnitudes {Vm and Rm) which were obtained with 
our broad band filters do not correspond to a standard V 
or Rkc (also K) exactly. Finally, as mentioned in Section 
4, there was an uncertainty of approximately ±0.015 mag 
dependent on the image position on each CCD chip. The re- 
sult from comparing the K magnitudes of Hughes & Wood 
(1990) and the average MO A magnitudes had a similar scat- 
ter to that due to the intrinsic variability of the stars. We 
also calculated the probable K magnitudes by extrapolating 
the observed light curves backwards to the time of observa- 
tion of Hughes & Wood (1990). The scatter found in this 
case was larger and most probably due to the imperfect pe- 
riodicity of Mira variables. Since the resulting relation was 
nearly the same as equation (14) that expression was used 
in the following study. 

Now we estimated the coincidence of our Km conver- 
sion with formula (14). In Fig. 12, we show a {{Vm) — 
{Rm), {Rm) — Kshv) plot for 19 oxygeu-rich stars (open cir- 
cles) and 11 carbon-rich stars (filled circles) which were com- 
mon in Hughes & Wood (1990) and our catalogue, and show 
the well-determined periods by our analysis. 19 oxygen-rich 
stars were used to derive formula (14) which is indicated by 
the solid line. It is obvious that the conversion formula (14) 
fits the oxygen- rich Miras (O-Miras) {o ~ 0.15 mag) well 
but does not fit the carbon-rich stars (cr ~ 0.34 mag). The 
dashed line indicates another conversion formula using 11 
carbon-rich stars: 



{Rm) - Kshv = 3.16 * {{Vm) - {Rm)) + 0.22 



(15) 



We also show the {[ogP,Km) diagram for these 19 
oxygen-rich stars and 11 carbon-rich stars in Fig. 13. The 
two lines in this figure were the sequences of O-Miras from 
previous studies (Feast et al. 1989; Hughes & Wood 1990). 
Our converted magnitudes Km for 19 oxygen-rich stars 
agreed well with these classical relations {o ~ 0.18 mag). 
In contrast, those for 11 carbon-rich stars appear below the 
lines (a ~ 0.25 mag). According to the results from previ- 
ous observations, the PL relation for O-Miras (oxygen-rich 
Miras) and C-Miras (carbon-rich Miras) are similar in the 
A'-band (Feast et al. 1989). In Fig. 12 and 13, we conclude 
that our conversion formula (14) makes a good enough fit 
to the oxygen-rich stars, but does not fit to the carbon-rich 
stars. 

We consider here why such a disagreement arises for 
the carbon-rich stars. As is well known, there are roughly 
two kinds of spectral types, M and C, for long-period vari- 
able stars depending on the relative abundance of oxygen to 
carbon in their atmospheres. This causes a remarkable dif- 
ference in colour tendency (Smak 1964). The B-band {Vm- 
magnitude in this paper) flux of oxygen-rich (M-type) stars 
was smaller than that for carbon-rich (C-type) stars. This 
was due to the existence of some strong TiO bands resulting 
from their rich abundance of oxygen. Therefore carbon-rich 
stars are bluer than oxygen-rich stars at the same effective 
temperature, which means that the value of {{Rm) — Kghv) 
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Figure 12. Comparison of {{Rm) — Kg^v) a^nd the MOA Colour 
((Vm) — {Rm,))- Kghv is the magnitude as tabulated in Hughes 
& Wood (1990). Open circles indicate the oxygen-rich stars and 
filled circles indicate the carbon-rich stars. The solid line indi- 
cates the Km conversion formula (14), while the dashed line in- 
dicates the least squares fitting for the 11 carbon-rich stars (for- 
mula (15)). The solid line makes a good fit to the oxygen-rich 
stars, while the 11 carbon-rich stars fall on the another steeper 
(dashed) line. 
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Figure 13. /Cj^^-magnitude versus Period for 19 oxygen-rich 
stars and 11 carbon-rich stars which were identified with Hughes 
& Wood (1990) as the same stars. The lines are the O-Mira se- 
quences obtained by previous studies (Feast et al. 1989; Hughes 
& Wood 1990). 



for carbon-rich stars is greater than that for oxygen-rich 
stars even if the colour index {{Vm) — {Rm)) is the same (see 
Fig. 12). This explains why the JC-magnitudes for carbon- 
rich stars were less luminous in Fig. 13, when we derived the 
conversion scale. Km, using only oxygen-rich stars. 

The above results were derived from the data by Hughes 
& Wood (1990). Within these stars (21 O-Miras and 13 
C-Miras), the if-magnitudes of nlmc3-306907 and nlmc3- 
288533 were reported in Table 4 of Trams et al. (1999). The 
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Table 2. /i'-magnitudes reported by other observations 



Star JD _fC/mag sources 

nlmc3-306907 2446341 10.96 Hughes and Wood (1990) 

2450788 10.35 Trams et al (1999) 

nlmc3-288533 2446341 10.58 Hughes and Wood (1990) 

2450782 10.65 Trams et al (1999) 



remarks 

phase fa 0.8 

near a light maximum 

phase ~ 0.8 

near a light maximum 



former is an oxygen-rich star and the latter is a carbon-rich 
star. Their results were shown in Table 2 together with Kshv 

Peak-to-peak modulation was found in the light curve 
of nlmc3-288533. The magnitudes tabulated in Table 2 were 
brighter than {Km), 11.35 and 11.20, of each star. The dif- 
ferences of magnitudes at the light maxima and {Km) were 
1.00 and 0.55 mag respectively. Such differences can be ac- 
cepted as the light variations of this type of variable stars. 

When we adopt the /f-magnitude by Trams et 
al.(1999), Kt, in Fig. 12, the position of nlmc3-306907 
moves upward and seems to belong to the carbon-rich stars. 
However, this star should not be regarded as a typical 
oxygen-rich star because of the unusual colour variation it 
showed as mentioned in Section 9.4. The position of nlmc3- 
288533 in Fig. 12 will be slightly lower than the current 
position. 

Figure 14 shows a similar diagram using the K- 
magnitude (Ks : ~ 2.15A^m) of th e DENIS catalogue 



(tittp://cdsweb.u-strasbg.fr/dcnis.html). Among 146 sam- 
ples, 139 stars have found, and JC-magnitudes were given 
for 119 stars in the DENIS catalogue. The correlation be- 
tween Ks and Km is given by the equation. 



Ks = 0.9(±0.09) *Km + 0.73(±1.07) 



(16) 



with a standard deviation of 0.47. Such a scatter may be ex- 
plained as the results of intrinsic light variability. In Fig. 14, 
oxygen-rich and carbon-rich stars were located in the lower 
and the upper regions as seen in Fig. 12. The star nlmc3- 
306907 was again located in the carbon-rich star region. No 
Ks data for nlmc3-288533 has been found. If nlmc3-306907 
is a typical oxygen-rich star, this indicates that some oxygen- 
rich stars may exist at the upper side of the oxygen-rich re- 
lation in Fig. 12. Therefore brighter oxygen-rich stars than 
given by equation (14) may exist. 

We show a (log P, Km) diagram for our sample of 146 
stars in Fig. 15. In this figure, the solid diagonal line is the 
LMC PL relation for O-Miras of Hughes & Wood (1990) and 
the dotted line shows the track for O-Miras of Feast et al. 
(1989). Three groups, designated A, B and C are evident 
in Fig. 15. They are discussed in detail in Section 11. The 
stars in the most populous group lie close to these lines. 
nlmc3-306907 is located in the same region as carbon-rich 
stars. Another star, nlmc2-64521, is also located below the 
established relation. Therefore it is clearly hard to decide 
the all star below these lines are carbon-rich stars. 

One of the 13 carbon-rich stars, nlmc2-64540, lies on 
these lines as an exception, and the remaining 12 carbon-rich 
stars fall below these lines. The slope of -ffs/iu-magnitude ver- 
sus MOA Colour (see Fig. 12 and formula (15)) for carbon- 
rich stars is steeper than that for oxygen-rich stars, so these 
lines intersect. This wEis why this exceptional carbon-rich 
star was located near the line of formula (14) and lies on the 
sequences for O-Miras in Fig. 15 (also in Fig. 13). 

We note it was possible that at least 50 ~ 60 carbon-rich 
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Figure 14. The {{Vm) ~ (Rm), (Rm) - Ks) relation for identi- 
fied 119 stars. Ks indicates the i^-magnitude from the DENIS 
catalogue (see text). 



stars were included in the sample of 146 stars because the 
ratio of oxygen-rich (the sample of 21 stars) and carbon-rich 
(the sample of 13 stars) stars was approximately ~ 62% and 
~ 38% respectively. These could generate the extra scatter 
in the (log P, Km) plot. 

When we converted the MOA data into K magnitudes 
by using equation (15), most of stars including 11 carbon- 
rich stars (around group C in Fig. 15) were found near the 
empirical LMC PL relation for C-Miras, but some of them 
remained in the fainter domain of the diagram. The isolation 
of the short period variable stars (seven -I- mark stars and 
a © star, defined in the next section) became clearer in this 
case. 



11 DISCUSSION 
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Figure 15. The (log P, Km) diagram for 146 long-period red variable stars. Two lines are the (logP, ii") relations for O-Mira sequences 
by the previous studies. The stars in the most populous group of our samples lie on the O-Mira sequence. One exceptional carbon-rich 
star also lies on the O-Mira sequence, while most of carbon-rich stars fall below these sequences. The + marks (seven samples) indicate 
the noteworthy stars with short period, and the ® indicates nlmc3-214266 (see text). An explanation of groups A, B and C is given in 
the text. 



11.1 log P — Km diagram 

It is difficult to estimate if-magnitudes from the large piio- 
tometry database as it was originally designed to detect mi- 
crolensing events. However we have now succeeded in esti- 
mating A"-magnitudes using the properties of oxygen-rich 
and carbon-rich stars as described in the previous Section. 
Based on these results, we can investigate the (log P, Km) 
diagram more carefully. 

As described in Section 10, our conversion of (Vm) and 
(Rm) into isTm-magnitude could not be applied for carbon- 
rich stars and they are located below the main sequence of 
Figures 13 and 15. Therefore the scatter in the (log P, Km) 
plane (Fig. 15) must be associated with these stars. How- 
ever, even taking the above fact into account, the scatter in 
the plot is too large to be described by one sequence. To 
demonstrate this point, in Fig. 16 we show a histogram of 
vertical distance from the line of O-Mira sequence by Hughes 
& Wood (1990): 

K = -3.86 * (log P - 2.4) + 11.5 (17) 

The error bars are equal to the square root of the number 
of stars in each of the three groups indentified in Figure 16. 

A vertical separation of zero means it was on the Hughes 
& Wood (1990) track in Fig. 15. The most populous group 
,A, lies at this position. This group contained the stars with 
largest amplitude. There was a small peak (group B, eight 



stars) around -1-1.5 mag (corresponding to a leftward shift 
in period in Fig. 15). The group B stars are relatively small 
amplitude variables. Another peak which is relatively broad 
can be identified at around —0.6 mag (group C, correspond- 
ing to a rightward shift in period in Fig. 15). We define 
each group in Fig. 16 using the following limits; the stars 
distributed between and 0.6 mag as group A, the stars 
distributed greater than 0.9 mag are as group B, and stars 
distributed between —0.9 and —0.3 mag as group C. 

In Fig. 17 we show the amplitude distribution, SRm for 
the stars in each group. In Fig. 18 we show the (log P, SRm) 
diagram with the same symbols as Fig. 15. In Fig. 17, it 
is clear that most of largest amplitude {5Rm > 2.2 mag) 
variable stars are group A stars, which provides strong proof 
that these stars are Miras. 

In contrast, the amplitudes of stars in group B were 
smaller {SRm < 1-4 mag) than most of the others as shown 
in Fig. 17 and Fig. 18. This group includes the seven short 
period stars and nlmc3-214266. This group is located on the 
very short period side (log P < 2.1) of the (Km, log P) plane 
(Fig. 15) . The group of seven short period stars is located in 
similar positions in the period-colour or period-magnitude 
diagrams (Fig. 7, 8, 9). All of these stars showed relatively 
inferior periodicity: the mean value of 9 was ~ 0.43, com- 
pared with the mean of 6 ~ 0.38 for the stars in group A. 
We obtained periodic variable stars via our selection criteria 
(Section 5), but we did not strictly eliminate stars whose pe- 
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Figure 16. The histogram of the vertical distance of i^-magnitude from the O-Mira sequence of Hughes & Wood (1990) in Fig. 15. The 
146 star sample is used. The highest peak at 5K = contains stars agreeing with the Mira (log P,i^) relation by Hughes &: Wood (1990). 
We classified these stars as group A. The stars separated by +1.5 mag from group A were defined as group B. The stars separated by 
—0.6 mag were defined as group C (see text). The error bars beside each group are described in the text. 
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Figure 17. The distributions of pulsation amplitude in Rm for each group. The upper panel shows stars on the most populous sequence 
(i.e., group A). The middle panel shows the distribution for the stars in group B which were located on the shorter period side of group 
A. The bottom panel shows that for the stars in the group C sequence which were located on the longer period side of group A. 



riodicities were relatively poor. So the stars in group B were 
somewhat semi-periodic compared to the stars in group A 
(see Fig. 5). 



11.2 The Group B Stars 

To examine K magnitudes of the short period variable stars, 
we compared our estimates with JHK observations made by 
the 1.9 m telescope at the South African Astronomical Ob- 
servatory. Star identification, magnitudes and date of obser- 
vations are listed in the first five columns of Table 3. Mean 



K magnitudes derived by equation (14) are tabulated in the 
6th column. The first three stars of the table are stars (a), 
(b), and (g) of group B, and the last one was the proba- 
bly non-AGB star (® in Fig. 15) of group B (see Table 1). 
The differences in K mag obtained at the SAAO and de- 
rived averaged Tf-magnitude, (Km), were 0.33, greater than 
1.02, 0.06, and 0.04, respectively. Except for the second star, 
nlmcl-239264, the differences do not seem remarkable. We 
predicted successive light maxima of these four stars using 
the observed light curves, and have found that the SAAO 
observations were done near a predicted light maximum for 
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Figure 18. The period-ampHtude {SRm) diagram. The symbols are the same as Fig. 15. Group B stars have been found in the low 
amplitude domain. 



the star nlmcl-324832, and near their predicted hght minima 
for the other three stars. If their hght variations were suf- 
ficiently regular, the SAAO if-magnitudes of nlmcl-324832 
had to be brighter than 10.47. This contradicts the observed 
brightness, so that the regularity of the light variation should 
be suspected. For the other three stars, /f-magnitudes were 
fainter than their {Km} - The observed magnitude for nlmcl- 
239264 was fainter than our adopted mean magnitude, and 
the difference was 1.02 mag. Such a difference looked too 
large when compared with the other two stars. However, 
the true difference between Ksaao and the mean K would 
be less than 1.02 mag because the adopted mean magnitude 
should be brighter than true mean magnitude. Therefore 
the difference may be acceptable with the true amplitude 
in Rm greater than 1.38 mag. The question remains as to 
whether these differences originate from the uncertainty in 
(Km) or the semi-regular nature of these variables. As indi- 
cated in Fig. 14, the Ks magnitudes for seven stars in group 
B are were located along the line of the 0-Miras. This tends 
to support the use of equation (14) in estimating the K- 
magnitudes. Further analysis of our database will resolve 
this question. 

Amongst the eight stars of group B, nlmc3-214266 
should be a massive star, not a star evolved from the red 
giant branch or the first red giant branch as pointed out in 
Section 9.1. We examine the properties of the other seven 
stars. The differences between the brightness of the group 
B stars and of the groupj4 stars is evident in Fig. 15. The 
faintest star, nlmcl-1189Q4, was equally bright in Ks with 
the other six stars of group B. The amplitude of nlmcl- 
239264 was larger than the other six stars, but the position 
on the period-luminosity diagram was not peculiar. We have 
to study the nature of this star by comparison with more 



abundant samples to determine whether or not this star is 
different from the others. 

In Fig. 15, the stars in group B were located parallel to 
the sequence of group A (corresponding to the Mira sequence 
of Hughes & Wood (1990)) and separated by (51ogP ~ 0.4. 
The seven stars in group B were located at a similar posi- 
tion to the short-period red (SP-red) stars which were re- 
ported by Feast & Whitelock (2000) analyzing Galactic Mi- 
ras. While the SP-red stars had large amplitude and Mira- 
like properties, our group B stars had small amplitudes ex- 
cept nlmcl-239264 and showed less regular variations, differ- 
ent from the Galactic Miras. The definition of the Galactic 
Miras is that their amplitude in standard 1/-light should be 
larger than 2.5 mag. It is a natural classification because 
there are few stars with such amplitudes. However, such a 
gap was not distinct in /-band observations of long period 
variables in the LMC (Whitelock, 1997). If we try to define 
Miras in the LMC by Fig. 17, the boundary might be at 
around 2.2 mag, though the small sample size makes that 
inconclusive. Therefore we discuss the long period variables 
of the LMC with no sharp distinction between Miras and 
semi- regulars. 

The theoretical models of long-period variables predict 
the ratio of fundamental to first or second overtone period 
should be S log P ~ 0.3 — 0.4. The stars in group B were very 
similar to those compared with the higher mode pulsators by 
Wood and Sebo (1996). Higher mode pulsators are also ex- 
pected to have smaller amplitudes than fundamental mode 
pulsators. The smallness of the amplitudes of these stars 
in Fig. 17 and Fig. 18 supported this point. Therefore if the 
stars in group A were fundamental pulsators, the seven stars 
in group B could be first or second overtone pulsators. On 
the other hand, as described above, these stars might be the 
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Table 3. JHK magnitude by SAAO 



Star 


J 


H 


K 


JD-2450000 




6Rm 


nlmcl-324832 


12.18 


11.12 


10.80 


2037.3 


10.47 


0.68 


nlmcl-239264 




12.25 


11.83 


2037.3 


10.81 


1.38 


nlmc2-19377 


12.43 


11.36 


10.96 


2037.3 


10.90 


0.73 


nlmc3-214266 


11.74 


10.92 


10.75 


2037.3 


10.71 


0.71 



Note: J, H, K magnitudes are measured to within ± 0.03 mag. 



counterparts of short-period red (SP-red) stars reported by 
Whitelock, Marang, and Feast (2000). These authors point 
out that the kinematics of such stars were different from the 
blue stars with similar periods. Such differences indicate a 
considerably different stellar age and imply a difference of 
evolutionary process due to a distinct chemical composition. 
A further development in the theory of stellar structure is 
needed before it will be possible to judge whether these seven 
stars are counterparts of SP-red stars (and are higher mode 
pulsators) or not. 

To study these stars further we plotted, in Fig. 19, the 
{Vm) — (Rm) colour agaiust K,n magnitude with the symbols 
being the same as in Fig. 15. The seven stars of group B are 
located at the reddest and highest luminosity region where 
longer-period variables are also located. This suggests that 
the radius of these stars is very similar to that of longer- 
period stars. The difference in the period of stars of similar 
radius is explained by the difference in the mass Af and the 
pulsation property Q. The estimate of the pulsation prop- 
erty for these stars has not been completed because of the 
uncertainty in the pulsation theory, but the effect of the 
chemical composition and the stellar mass on Q looks neg- 
ligible for a given mode (Barthes and Tuchman 1994). The 
difference in the period-temperature relation found in the 
globular clusters (Whitelock 1986) has been resolved using 
different chemical compositions, but the difference in the pe- 
riod of the stars of which the colour is the same cannot be 
explained this way. The difference in the stellar mass and 
the mode of pulsation should be studied to explain the dif- 
ference in the period. When the distribution of periods with 
stellar physical properties is continuous, it is reasonable to 
suppose that this is due to changes in stellar mass or other 
stellar parameters such as chemical composition. Since the 
periods of these seven stars are separated from their long- 
period counterparts, the enhancement of a different mode 
is suggested. Further progress in data reduction is required 
to examine the existence of the gap in the periods of this 
luminous red domain of the colour- A" diagram. 

When we consider the excitation of different modes, we 
have to study the reason for such a mode selection. In the 
classical Cepheids and RR Lyrae stars, the higher mode 
pulsators have been found on the blue side, or small stel- 
lar radii side, of the colour-magnitude diagram. Therefore it 
seems difficult to compare directly the group B stars with 
higher mode pulsators of the Cepheid Instability Strip. On 
the other hand, it may be interesting to compare the sur- 
face gravity between them. The surface gravity of the higher 
mode pulsators is high for stars in the Cepheid Instability 
Strip. This implies that the surface gravity of the group B 
stars might be also high. The critical frequency of the stel- 
lar atmosphere (i.e., the boundary frequency of reflection or 
transmission at the stellar surface) must be high for a high 
surface gravity. The enhancement of high frequency, or short 



period pulsation in group B stars suggests high stellar mass 
when we suppose the same situation for the transmission 
and reflection for all pulsating stars. 

11.3 The Group C Stars 

There are 46 stars in group C. Many carbon-rich stars might 
be included in this group. 12 group C stars are the same as 
the carbon-rich stars tabulated in Hughes & Wood (1990). 
There was a possibility that about 50 ~ 60 stars among the 
sample of 146 are carbon-rich stars as mentioned in Section 
10. These carbon-rich stars were estimated to be too faint 
in our Km conversion. Also, in the amplitude distribution 
(Fig. 17), there was no obvious difference between group A 
and group C. We conclude most group C stars were carbon- 
rich and should be included in group A. However, we need 
to distinguish oxygen-rich stars from carbon-rich stars in 
order to investigate whether a long-period sequence exists, 
and if the strip structure is real or not. Future study of the 
relationship using a B-band filter should help resolve this 
issue. 

11.4 Strip Structure 

Wood (2000) presented the sequential structure of the PL 
diagram using a large sample. Because our sample size is 
small, it was hard to confirm this structure, except for the 
existence of a separate group of red short-period stars. How- 
ever, it should be noted that these red short-period stars 
were of small amplitude {5Rm) and less regular (small 6^). 
These will be useful characteristics in the study of the nature 
of these stars. In the period-iC diagram, the stars, nlmc3- 
214266 and nlmcl-239264, were also contained in group B 
implying that it is necessary to investigate the properties of 
individual star even in large scale photometry in order to 
understand the period-luminosity relation. The individual 
star discovered by Wood (2000) should also be checked to 
determine whether it is to be included in group B and the 
in the set of Galactic red SP stars. 



12 SUMMARY 

We selected long-period red LMC variables from the three 
year MOA database, and executed a careful period analysis 
using the PDMM and PDM codes. We selected 146 variable 
stars using criteria which focus on long-period red variables 
(Miras and SRs), and investigated their properties. Their 
observational data were tabulated. 

In the (logP, (Vm)) and (logP, (Rm)) diagrams, the lu- 
minosity of stars in the major cluster tended to decrease as 
the period increases. In the (log P,{Vm) — (Rm)) diagram, 
the colours were redder as the periods were longer. 
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Figure 19. The {{Vm) — {Rm), Km) diagram for 146 samples. The symbols are the same as Fig. 15. The seven stars with short period 
(+) are located at the reddest and highest luminosity region. 



We transformed Vm , Rm magnitudes to Km values using 
19 oxygen-rich stars, and presented the period-luminosity 
(PL) relation in the (logP, i^m) plane. Many of the stars in 
our sample, including the large amplitude stars, fell near the 
position corresponding to the O-Mira sequence by Feast et 
al. (1989) and Hughes & Wood (1990). A smaU number of 
stars (8) formed a separate group in the period-magnitude or 
period-colour diagrams. This group (group B) was found in 
the bright and short period region of the (log P, Km) plane 
compared with the O-Mira sequence. These stars showed less 
periodic variation, and six of them had small amplitudes. 

The separation in log P between them and the main 
group was ~ 0.4. They were on the red side of the period- 
colour diagram and were located at the reddest and bright- 
est tip of the colouT-Km diagram similar to the long-period 
members. If our conversion scale is adequate, it indicates 
that some of the stars pulsate in a short period compared 
with the stars which may have a similar radius. 

These facts were possibly explained by the assumption 
proposed by Wood & Sebo (1996): that those six stars with 
short periods and small amplitudes were higher mode pul- 
sators than the stars on the O-Mira sequence. Two partic- 
ular stars, nlmc3-214266 and nlmcl-239264, were included 
in group B. These stars showed different features from the 
other six variable stars. Their pulsation property possibly 
differs from the other six stars. 

In the (log P, Km) diagram, below the O-Mira sequence, 
many stars were located in a wide region. Since almost all 
identified C-Miras were found in this domain, we supposed 
that a large number of such variable stars were C-Miras 
which were more luminous than the magnitude indicated 
here. In order to reduce the systematic scattering in the 



(log P, Km) diagram, it was necessary to distinguish carbon- 
rich stars from oxygen-rich stars in advance. 

The difficulty in distinguishing oxygen-rich and carbon- 
rich stars comes from the fact that our blue filter covers 
both the standard B and V bands. It is expected that in the 
future these stars will be observed in the B or V bands which 
will distinguish between carbon-rich stars and oxygen-rich 
stars. Alternatively, using any passband which avoids the 
effect of TiO bands to estimate the bolometric magnitude 
of the stars could be used. Improvement in the magnitude 
calibration and confidence in the interpretation of the plural 
sequences is expected to follow from the analysis of further 
large observational databases. 
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